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Abstract 

We study long wavelength perturbations of neutral black p-branes in asymptotically flat 
space and show that, as anticipated in the blackfold approach, solutions of the relativistic 
hydrodynamic equations for an effective p + 1-dimensional fluid yield solutions to the vacuum 
Einstein equations in a derivative expansion. Going beyond the perfect fluid approximation, 
we compute the effective shear and bulk viscosities of the black brane. The values we obtain 
saturate generic bounds. Sound waves in the effective fluid are unstable, and have been 
previously related to the Gregory-Laflamme instability of black p-branes. By including the 
damping effect of the viscosity in the unstable sound waves, we obtain a remarkably good 
and simple approximation to the dispersion relation of the Gregory-Laflamme modes, whose 
accuracy increases with the number of transverse dimensions. We propose an exact limiting 
form as the number of dimensions tends to infinity. 



1 Introduction and Summary 



Black holes exhibit thermodynamic behavior, so it is natural to expect that their long wave- 
length fluctuations, relative to a suitable length scale, can be described using an effective 
hydrodynamic theory. Over the years there have appeared several different realizations of 
this idea, which differ in the precise set of gravitational degrees of freedom that are captured 
hydrodynamically (e. 5., only those inside a (stretched) horizon as in [Ij, or the entire gravi- 
tational field up to a large distance from a black brane spacetime as in [21 13] ) or in the kind 
of asymptotics (Anti-deSitter [2] or flat [3]) of the black hole/brane geometry. 

In this paper we focus on the hydrodynamic formulation developed recently for higher- 
dimensional black holes, including asymptotically flat vacuum black holes and black branes 
[3]. In this approach the effective stress tensor of the 'black brane fluid' is the quasilocal stress 
tensor computed on a surface i3 in a region that is asymptotically flat in directions transverse 
to the bran^. The equations of stress-energy conservation describe both hydrodynamic 
(intrinsic) fluctuations along the worldvolume of the brane, and elastic (extrinsic) fluctuations 
of the brane worldvolume inside a 'target' spacetime that extends beyond B. Thus the 
dynamics of a black p-brane takes the form of the dynamics of a fluid that lives on a dynamical 
worldvolume. This is referred to as the blackfold approach. 

In this paper we only study the intrinsic, hydrodynamic aspects of the brane. The world- 
volume geometry, defined by the surface B at spatial infinity, is kept flat and fixed. Fluctua- 
tions of the worldvolume geometry are non-normalizable modes, so the extrinsic worldvolume 
dynamics decouples. With this simplification, the set up is very similar to the fiuid/AdS- 
gravity correspondence of p], which we follow in many respects. The main difference is that 
we consider vacuum black brane solutions, with no cosmological constant and with different 
asymptotics. 

The quasilocal stress tensor of a neutral vacuum black brane, with geometry equal to the 
n -|- 3-dimensional Schwarzschild-Tangherlini solution times M^, is that of a perfect fiuid with 
energy density p and pressure P related by the equation of state 



brane could also be boosted and thus have a non-zero velocity field along its worldvolume. 
In a stationary equilibrium state, the temperature and the velocity are uniform. We study 
fiuctuations away from this state where these quantities vary slowly over the worldvolume. 
Their wavelength is measured relative to the thermal length T~^, so for a fiuctuation with 
wavenumber k the small expansion parameter is 



Since for a vacuum black brane the temperature is inversely proportional to the thickness of 
the brane, rg, this can be equivalently expressed as krQ <^ 1. 

'^In the following, asymptotic flatness always refers to directions transverse to the brane. 
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To leading order in this expansion we obtain the hydrodynamics of an effective perfect 
fluid, which refs. [H EJ [5] have used to derive non-trivial results for higher-dimensional black 
holes. At the next order the stress tensor includes dissipative terms. For the purely intrinsic 
dynamics, these are the shear and bulk viscosities, rj and C. In contrast to [2j, our fluid is 
not conformally invariant so C 7^ is expected. 

By analyzing long wavelength perturbations of the black brane and their effect on the 
stress tensor measured near spatial infinity we obtain 

where s is the entropy density of the fluid, z.e.,l/4G times the area density of the black 
brane, and 

= 1 a.4) 

" dp n + 1 ^ ^ 

is the speed of sound, squared. 

Written in the form ()1.3p , these values for r] and C saturate the bounds proposed in ^ and 
[7]. The result for the shear viscosity is not too surprising: rj can be argued to depend only 
on the geometry near the horizon and its ratio to s is universal for theories of two-derivative 
Einstein gravity [H |8| (see also [9]). The bulk viscosity, instead, does depend strongly on 
the radial proflle transverse to the bran^ so the saturation of the bound is presumably less 
expected. Note, however, that these black branes have different asymptotics than in all the 
previous instances where the effective viscosities of black branes have been considered. In 
particular, these black branes presumably are not dual to the plasma of any (local) quantum 
field theory. In any case it is worth emphasizing that our computations are for the theory 
with the simplest gravitational dynamics: R^u = 0. 

The imaginary speed of sound (jl.4p implies that sound waves along the effective black 
brane fluid are unstable: under a density perturbation the fluid evolves to become more 
and more inhomogenous. Since this means that the black brane horizon itself becomes in- 
homogeneous, ref. [3] related this effect to the Gregory-Laflamme (GL) instability of black 
branes [luj^ . Then (|1.4|) implies a simple form for the dispersion relation of the GL unstable 
modes uj{k) = —iQ{k) at long wavelength: J7 = k/^/n + 1 + 0{k^), i.e., the slope of the 
curve ^}{k) near A: = is exactly (and very simply) determined in the unstable-perfect-fluid 
approximation. 

Using our results for ij and C we can include the viscous damping of sound waves in the 
effective black brane fluid. The dispersion relation of unstable modes becomes 

^ = y^(^-^^kro) , (1.5) 
\/n + 1 V n\/n +1 / 

which is valid up to corrections oc A;'^. Figure [J compares this dispersion relation to the 
numerical results obtained from linearized perturbations of a black p-brane. Zooming in on 



^For instance, in the membrane paradigm the bulk viscosity on the stretched horizon for a generic black 
hole turns out to be negative. Our result (|1.3[) is instead positive. 

^This connection had also been made for black branes with gauge theory duals in 
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Figure 1: Left: dispersion relation eq. (jl.Sp . for unstable sound waves in the effective 

black brane fluid (normalized relative to the thickness rg). Right: for the unstable 

Gregory-Laflamme mode for black branes (numerical data courtesy of P. Figueras). For 
black p-branes in D spacetime dimensions, the curves depend only on n = D — p — 3. 



small values of Zero, the match is excellent. When kr^ is of order one we have no right to expect 
agreement, but the overall qualitative resemblance of the curves is nevertheless striking. The 
quantitative agreement improves with increasing n and indeed, as figure [2] shows, at large 
n it becomes impressively good over all wavelengths: for n = 100 the numerical values are 
reproduced to better than 1% accuracy up to the maximum value of k. Although the extent of 
this agreement is surprising, we will provide some arguments for why the fiuid approximation 
appears to be so successful as n grows. 

Thus, the effective viscous fluid seems to capture in a simple manner some of the most 
characteristic features of black brane dynamics. We believe this is a significant simplification 
from the complexity of the full Einstein equations. 

The outline of the rest of the paper is as follows: the next section contains the bulk of the 
calculations of the paper for a generic hydrodynamic-type perturbation of the black brane. 
We highlight the differences with the analysis of [2], in particular at asymptotic infinity, and 
compute the values (jl.3p for the effective rj and Q. Section [3] relates the linearized damped 
sound-mode perturbations of the fluid to the Gregory-Laflamme perturbations of the black 
brane. We examine the conditions that can lead to the surprising quantitative agreement of 
the dispersion relation at large n, and we propose its exact form as n — )• oo. We close in 
section [J] with an examination of the differences with other fluid-like approaches to the GL 
instability, and a discussion of our results within the context of the blackfold approach. 



2 Hydrodynamic perturbations of black branes 

In this section we study general perturbations of a vacuum black p-brane with slow variation 
along the worlvolume directions of the brane. Up to gauge transformations, they are fully 
determined by the boundary conditions of horizon regularity and asymptotic flatness at spa- 
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Figure 2: Dispersion relation 0,{k) of unstable modes for n = 100: the solid line is our 
analytic approximation eq. (jl.Sp ; the dots are the numerical solution of the Gregory-Laflamme 
perturbations of black branes (numerical data courtesy of P. Figueras). 



tial infinity. Most of our analysis is very close to the study of hydrodynamic perturbations of 
AdS black branes, but there is an additional complication in the study of the perturbations 
at asymptotic infinity. Nevertheless, we are able to find the complete explicit form of the 
perturbed solution for a generic hydrodynamic flow to first order in the derivative expansion. 

Readers who do not need or want the technical details of the calculation of the perturbed 
solution and the viscous stress tensor can safely skip to section [3l 



2.1 Preliminaries 

2.1.1 Black branes and their effective stress tensor 

The black ]5-brane solution of vacuum gravity mD=p + n + 3 dimensions is 

Vab + :^UaUbj da" da'' + ^ + r'^d^l^^ , (2.1) 

with a = 0, 1, . . . The solution is characterized by the horizon radius rg (or brane 'thick- 
ness') and the worldvolume velocity with u°'v}'riab = —1. It is asymptotically flat in the 
directions transverse to the worldvolume coordinates a"". We can associate to it a stress- 
energy tensor measured at spatial infinity. There are several possible definitions of this stress 
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tensor that would be equivalent for calculational purposes, but for conceptual reasons the 
most convenient for us is the quasilocal one of Brown and York [12j. We consider a boundary 
surface at large constant r, with induced metric /i^jy and compute 

T^r^ = ^ {k,, - h,uK - {Kj^J - V^^°^)) , (2.2) 

where K^^, is the extrinsic curvature of the surface and we perform a background substraction 
from flat spacetime. 

The geometry of the boundary surface for ()2.ip is R^'P X 

"YYg ^j^j introduce per- 
turbations with wavelengths much longer than the size rg of the S^~^^ at the horizon. The 
deformations of this sphere all have large masses ~ 1/ro and therefore decouple. Thus the 
SO{n + 2) symmetry of S""^^ is preserved and, in an appropriate gauge, the metric will re- 
main a direct product with a factor of this sphere. We integrate over the sphere to obtain 
the stress tensor for the black p-brane 

Ta, = [ T^r^ . (2.3) 

We regard this stress tensor as living on the worldvolume of the brane, i.e., the p + 1 
extended directions of the boundary. The worldvolume metric results from the asymptotic 
form of the boundary metric, which in our case is the Minkowski metric 

hab r]ab . (2.4) 

A main advantage of using the quasilocal stress tensor is that the Gauss-Codacci equations 
for the constant-r cylinder imply d"'Tab oc Wb^ so imposing the Einstein equations in vacuum 
it follows that the stress tensor is conserved 

d'^Tab = . (2.5) 

The stress tensor for the solution (|2.ip has the perfect fluid form 

Tab = PUaUb + PPab , Pab = Vab + UaUb (2.6) 

with energy density and pressure 

p=-(n+l)P=(n+l)^^. (2.7) 
The horizon area allows to associate a local entropy density to this effective fluid 

and all the thermodynamic functions can be expressed as functions of the temperature 

„ n , ^ 

T = . 2.9 

47rro 
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We can equivalently use T or tq as the variable that determines local equilibrium. In this 
section we will mostly use tq for notational simplicity. 

We will be interested in preserving regularity at the horizon. This is manifest if instead 
of the Schwarzschild coordinates in (j2.ip we use Eddington-Finkelstein (EF) ones, 

a^^a'^- u'^r. , = j ^3^^^^ ' (^-^O) 

such that 

ds^ = - ^1 - UaUbda'^da^ - 2u„dcr"dr + Pabda^da^ + r'^dnl^i . (2.11) 
2.1.2 Perturbations 

We promote the thickness and velocity parameters to collective fields over the worldvolume, 
so 

dsfo) = - ( 1 - ) ua{a)ub{a)da^da^ - 2ua{a)da''dr + (r/ab + Ua{a)ub{a)) da^da^ 



+r^dnl^, , (2.12) 



where ro{a) and u'^{a) are assumed to vary slowly relative to the scale set by tq. In this 
paper we expand them to first order in derivatives, which we keep track of through a formal 
derivative-counting parameter e. With non-uniform and u"", the metric (I2.12p is not Ricci 
flat so we add to it a component with radial dependence 

ds^ = dsf^-^ + ef^,y{r)dx^'dx'' + 0{e^) . (2.13) 

We choose a gauge in which dr is a null vector with normalization fixed by the radius r of 
S'^+i, so that 

/rr = 0, /q^ = 0. (2.14) 

With this choice the sphere ^"""^^ can be integrated out. 

Demanding that (I2.13P satisfies the vacuum Einstein equations to first order in e results 
into a set of ODEs for f^y{r). These will be solved subject to regularity at the horizon 
r = ro, which is easily imposed as a condition of metric finiteness in EF coordinates, and to 
asymptotic flatness, to which we turn next. 

2.1.3 Asymptotic infinity 

The asymptotic behavior of our spacetimes introduces an important difference relative to the 
perturbations of AdS black branes. For the latter, the calculations can be performed in their 
entirety in EF coordinates in which dr is a null vector. By taking large values of r in these 
coordinates one approaches null infinity, but in AdS this is the same as spatial infinity. The 
AdS boundary is always a timelike surface. However, in our asymptotically flat space, null 
and spatial infinities differ. 
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We are ultimately interested in computing the quasilocal stress tensor on a timelike bound- 
ary of spacetime endowed with a non-degenerate metric. But if we approach null infinity, 
the boundary metric will be degenerate and it is unclear whether the quasilocal stress tensor 
is well defined there. Instead, it seems more appropriate (and is definitely unproblematic) 
to compute the stress tensor at spatial infinitjo. For this purpose EF coordinates are very 
awkward and it is much more convenient to switch back to Schwarzschild-like coordinates 
{r, t, cr*} at large r. 

Thus we will work with two sets of coordinates: EF ones, in which horizon regularity is 
manifest, and Schwarzschild coordinates, in which spatial infinity is naturally approached. 
We need to provide the change of coordinates that relates them, extending the inverse of 
(j2.10p to include 0(e) terms. The correction is naturally guessed by recalling that and rg, 
which appear in the transformation ()2.10p . now depend on the EF coordinates. Thus, 



a-^a'' + u%v, a') / f'' , (2.15) 

' 1 — [ro{v, a^)/r)"- 



or more explicitly, 



ro V l-(ro/r)"y 

a' ^ a' + e {{t + u)dyu' + a^dju') + 0{e^) . (2.16) 

2.2 Solving the perturbation equations 

At each point we choose coordinates centered on that point and go to an (unperturbed) local 
rest frame. In EF coordinates the velocity perturbation is 

u^'ia) = l + 0{e^), u\a) = ea''dau'{0) + O{€'^). (2.17) 

Note that since local velocities are small the constraint = —1 is automatically satisfied 
to the order we need. The other collective variable of the effective black brane fluid is the 
temperature T, or equivalently the thickness ro, which we perturb as 

ro(fT) = ro(0) + ea'^daroiO) + 0{e^) . (2.18) 

In the following we understand all quantities as evaluated at cr" = and thus denote 
dau'iO) dau\ ro(0) ro etc. 



^Presumably the appropriate notion of spatial infinity here is not Penrose's (which is just a point) but 
more along the lines of [13], which naturally allows a dependence along the boundary directions. Although our 
spatial infinity is not exactly the same as in 113. since instead of a hyperboloid we work on a cylinder where 
R^''' and 5""''^ scale differently at infinity, this is not a problem for us since we are integrating over S"^^ . It 
would be interesting, especially with a view to holography, to further formalize this notion of spatial infinity. 
Related remarks concerning holography in asymptotically fiat spacetimes have been made in [14j . 



7 



The metric (|2.13p is now 

p 

ds^ =2dvdr - f{r)dv^ + ^ daf + r'^dQl^^ 

- 2ea^dauMdr + g ^^o"'^"^"^o _ 2jl^^^da'dv + ef^,{r)dx''dx'' , (2.19) 



where we denote 

fir) = 1 - i . (2.20) 



The Einstein equations with a radial index, R^a = do not involve second derivatives 

and are constraint equations. Indeed they only involve the hydrodynamic fields ro and 
and not ff,,^, 

(n + l)dyro = -rod^i , diV^ = r^d^m, (2-21) 

so they are to be regarded as the equations of fluid dynamics, consistently with (|2.5p . We 
also verify this interpretation later. 

The remaining Einstein's equations are dynamical and we solve them to find f^y. The 
equations Rij = give 

dr {r'^^'fn/) = -2(n + l)r"%n,) , (2.22) 
which, requiring finiteness at the horizon, are solved by 

fijir) = Cij - 2d(^iUj) (^r^ - ^ log . (2.23) 

The integration constants Cij will be fixed later demanding asymptotic flatness. The equations 

dr = - (n + l)r^d,u, , (2.24) 

are solved by 

fvi = cf> + ^-d,Uir, (2.25) 

which are regular at the horizon for all values of the constants. Next, the equations from 
Rrr = and Rqq = are 

and 

dr (r"/™) = r'^d'u, + ^ (^E - - 2nr"-V.r , (2.27) 
which, assuming that eqs. (|2.21|) are satisfied, are solved by 

d fa 

jvr = Cot + 777 ; 7T3~ / ^ ' (2.28) 
2(n + 1 dr ^-^ r 

i 

and 

_ 2gW+(l-^g) ELi/n rj 

fvv — ; ^ 2ct,j. H ~Cnv ■ (2.29) 
n + 1 r" 
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Again these are regular at the horizon for all choices of the integration constants. Note that 
frj does not appear in Einstein's equations to first order in e and corresponds to a gauge 
mode. This, and the integration constants, will be fixed shortly. 

At this stage, for any hydrodynamic perturbation that solves the equations (|2.2ip . we have 
managed to construct a perturbed metric that is regular at the horizon. Next we must ensure 
that the solution remains asymptotically flat. Transforming to Schwarzchild-like coordinates 
using (j2.16p . we require that 

9ab = Vab + 0(r-") . (2.30) 

For the other metric components, we find that Qrr = 1 + 0(r~"), g^i = 0{r~"'), and gtr = 
0{r~^^^) when n > 1 [gtr = 0(logr/r) when n = 1), are enough to obtain a finite stress 
tensor. Recall also that all the metric components involving angular coordinates of S^~^^ are 
unaltered. 

Omitting details, we find that the conditions on gij and gtj fix 

=cg)=0. (2.31) 

In addition, the effect of c^j in gtj amounts to a global shift in the velocity field along the 
spatial directions of the brane, so in order to remain in a local rest frame we set 

ci'^=0. (2.32) 

Furthermore, if we perform the change 

t^t{l- tc„) , (2.33) 
then Cm — 2ct,r results in a global shift in the temperature, which we eliminate by choosing 

^VV 2C|;7- . (2.34) 

Asymptotic flatness in imposes a choice for c^- that singles out the slower fall-off of n = 1, 

c-vr = —QtTQ for n = 1 , c^r = for n > 1 (2.35) 

(note that the values of Strg and d^r^ at a" = are equal). 

Asymptotic flatness in these coordinates is a little delicate when n = 1 due to its slower 
fall-off, and to make it manifest we take an /^j gauge diverging at infinity. This is not 
necessary when n > 1 (and neither choice affects the calculation of the stress tensor). Thus 
we set 

r 

frj = —djTQ log — for n = 1 , frj=0 for n > 1 . (2.36) 
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Summarizing, we obtain 

g^j=Sij+ero^^^^\ogf, (2.37) 



n 

„a 



9tj = - ^^.^"dau, , (2.38) 



5n=e/.i(r) + 6^^, (2.40) 
5rr =/ + ef [-^^ + (log / - 2) J , (2.41) 

5« = -/ + e(^^^7;^ + ^*^olog/(^^--/J J , (2.42) 

((T° correspond to Schwarzschild coordinates here, so = t). This is the complete solution 
for the black brane metric that corresponds to a hydrodynamic perturbation that solves the 
equations ()2.2ip expanded around the origin of the local rest frame, a°- = 0. 

2.3 Viscous stress tensor 

We are now ready to compute the quasilocal stress tensor (j2.2p . The renormalization via 
background subtraction is simple and appropriate, since our metrics are infinitesimally close 
to the uniform black p-brane and their asymptotic boundaries can always be embedded in 
flat spacetime. Straightforward calculations give 



Ttt = + 1) (ro + ea'^daVoT , (2-43) 



25(jUj) - '^6ijd^Ui^ + 2 Q + — ^ ) 6ijd'^ui 



i'n+1'0 a a 



which are valid up to O(e^). One can easily check that the hydrodynamic equations dgT""^ = 
are indeed equivalent to the constraint equations (j2.2ip . 
Write now this stress tensor in the form 

Tab = pUaUb + PPab " C^Pab " 2r]aab + 0{d^) (2.44) 

where the expansion and shear of the velocity congruence are 

1 

The component Tu in (j2.43p determines the energy density, and requiring that the equa- 
tion of state (|l.ip holds locally uniquely identifies the pressure. Then we can write 

Tij = P5ij - eri (^29(,n,) - ^d^jd%^ - eC6^jd% (2.46) 



daU'^ , (Jab = Pa' ( 5(,nrf) - -P,d ) P\ • (2-45) 
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with 

Using (jl.4p and (|2.8p these can be rewritten as in (|1.3p . 

3 Damped unstable sound waves and the Gregory- Laflamme 
instability 

Our analysis in the previous section apphes to generic long-wavelength perturbations of arbi- 
trarily large amplitude. Let us now consider small perturbations of a static fluid of the form 

p^p + 5p, P^P + cl5p, u" = (1,0,...) ^ {l,6u'), (3.1) 
where Cg is the speed of sound, and with 

5p{t, a') = 5pe^^'+^^^''' , 5u\t, a') = 5u' e*^*+*'=^'^' . (3.2) 

We substitute these in the viscous fluid equations and linearize in the amplitudes 5p and iJu*, 
to find 

u5p + {p + P)ki5u' + 0{k^) = (3.3) 
ioj{p + P)5u^ + icl^dp + r]k^5u^ + /c^' ^ ^1 - r/ + hdu^ + 0{k^) = . (3.4) 

Applying our results above, any solution to these equations can be used to obtain an explicit 
black brane solution with a small, long- wavelength fluctuation of tq and n". If we eliminate 
6p we find that non-trivial sound waves require 

u-cl--t^i2(l-^\r^ + c\+ 0{k') = , (3.5) 



LU Ts \ \ p 

where k = \/kiki and we have used the Gibbs-Duhem relation p + P = Ts. This equation 
determines the dispersion relation u:{k). For a stable fluid with > 0, viscosity adds a 
small imaginary part to the frequency, which becomes complex and describes damped sound 
oscillations. Instead our effective fluid has imaginary sound-speed, eq. (II. 4p . so u is purely 
imaginary: sound waves are unstable. Writing 

u = —iO. (3.6) 

we solve (13.51) to flnd 



(3.7) 
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For the specific black p-brane fluid this yields the dispersion relation p.5|) . The connection 
between these unstable sound waves and the Gregory-Laflamme instability was pointed out at 
the perfect fluid level {i.e., linear in k) in [3J, and we have discussed it in the introduction 

Figures [1] and [2] show that our approximation (jl.Sp improves as n grows. In order to see 
how this might be justified, let us first rewrite the dispersion relation (jl.Sp in terms of the 
temperature T instead of rg, 



In principle, at any given n, both quantities tq and T define length scales that are paramet- 
rically equivalent. But if we vary n and allow it to take large values, then rg and T^^ ~ ro/n 
can differ greatly. We propose that in this case, T~^, and not rg, is the length scale that 
limits the validity of the fluid approximation, so the appropriate expansion variable for large 
n is k/T and not krQ. This may actually be natural since from the fluid point of view T has 
a clearer physical meaning than tq. In effect, we are proposing that when n ^ 1 it is more 
accurate to view the effective theory as describing very hot black branes, rather than very 
thin ones. 

The point of this exercise is that for large n the maximum values over which il. and k in 
(|3.8|) range are (A;/T)|niax ~ I/aA^ and {^/T)\ 

max ~ 1/?T.. So as Ti grows the frequency and 
wavenumber of unstable modes extend over a smaller range of k/T and fi/T. This strongly 
suggests that hydrodynamics can capture more accurately the dynamics of GL modes when 
the number of dimensions becomes very largely More precisely, if we write the corrections 
inside the brackets in (|3.8p in the form X]j>2 i aiid assume that the n-dependence 

of the coefficients aj is such that ajn^^^"^ — )• as n — t- oo, then the expansion in k/T, i.e., the 
hydrodynamic derivative expansion, becomes a better approximation over a larger portion of 
the curves Q{k). 

This is a relatively mild-looking assumption on the n-dependence of the higher-order 
coefficients in the expansion in fc/T|j and in particular is satisfied if the aj>2 remain finite as 
n — >• oo. But since we have not computed higher-derivative transport coefficients then, within 
our perturbative framework, we cannot prove its validity. However, since the numerical data 
appear to strongly support it, we conjecture that the truncation of the dispersion relation up 
to fc^-terms captures the complete dispersion relation at large n. More precisely, if we define 
a rescaled frequency and wavenumber. 



^Observe that the result (|1.5|) is independent of p. That this must be the case is clear from the outset in 
the GL analysis and also in our analysis of the Einstein equations. 

^This is similar in spirit, although not precisely equal, to the proposal in 15 that in the limit of large 
number of dimensions black holes are accurately described by fluid mechanics. 

Which, crucially, is not satisfied by the coefficient of the linear term inside the brackets in (|3.8p . 
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that remain finite as n — >• oo, then we propose that 



n = k\i 



AttT 



k 



) 



(3.10) 



is the exact limiting relation valid for all wavenumbers < k < AtiT. 

Note that the trmication of Vt{k) in (|3.8p appears to capture the zero- mode with = 
at a finite k = kcL- This is quite remarkable, since the viscous fluid equation (13.5P does 
not admit any zero-mode solution. The comparison with numerical data in figure [2] shows 
that the quantitative result for kcL-, although poor for small n, becomes excellent for large 
n. Further evidence for the validity of our proposal comes from the analytical value of the 
GL zero mode in the limit n — ?■ oo |16] 



This is the same as the limiting value for the zero- mode 'predicted' by (13. 9p . (I3.10p rl 

Presumably, by effecting the scaling (j3.9p in the full linearized perturbation equations of 
the GL problem one may prove (or possibly disprove) equation (j3.10p . 

4 Discussion 

Our analysis of the GL instability must not be confused with recent studies where a connection 
to the Rayleigh-Plateau instability of fluid tubes is made. In the latter approach, following 
a suggestion in |18j . refs. [151 fT9] related a d-dimensional black string in a Scherk-Schwarz 
compactification of Anti-deSitter space to a d — 2-dimensional fluid tube with a boundary 
with surface tension (see |20j). The Rayleigh-Plateau instability of the fluid tube arises from 
the competition between surface tension and bulk pressure. In contrast, our effective fluid 
does not have any boundaries so the instability is not of the Rayleigh-Plateau type, but 
rather one in the sound modes. Also note that our calculations in sec. [2] yield explicit black 
brane solutions to the Einstein equations (in vacuum) in a derivative expansion, something 
that, although expected to be possible in principle, at present cannot be realized for the fluid 
solutions in [15l [19] . 

We stress that our analysis is not a 'dual' solution of the GL instability problem: we have 
investigated the same perturbation problem as in [10] and explicitly solved it in closed analytic 
form in a derivative expansion. Since our approach does not require the perturbations to be 
small, it may even be used to study the non-linear evolution of the GL instability. 

One of our motivations has been to show explicitly how the effective theory of blackfolds 
of [3] can be systematically developed as a derivative expansion of the Einstein equations. 
Although we have done it only for the intrinsic aspects of blackfold dynamics, we have been 

*The relative difTerence between the results for kcL from the large-n subleading correction computed in 
[17) and from (II. 5p is equal to 1/n. This is precisely the size of the discrepancy observed in fig. (21 



kcL 



AttT 



(3.11) 
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able to: (i) derive in detail, starting from the 'microscopic' (full Einstein) theory, the lowest- 
order blackfold formalism that ref. [3] had developed following general principles; (ii) prove 
that the first corrections to the lowest-order formalism can be computed and result in pertur- 
bations of the black brane that preserve regularity of the horizon. The viscosity coefficients 
are determined precisely from this condition. 

In general, the worldvolume of a blackfold is dynamical and can be curved. Our calcula- 
tions in this paper can be regarded as being valid for fluid perturbations with a wavelength 
that, while longer than T~^, is much shorter than the typical curvature radius R of the black- 
fold worldvolume. In this case, the intrinsic and extrinsic dynamics decouple. Thus, for a 
curved blackfold our results for the GL instability are valid at most up to wavelengths smaller 
than R. At longer wavelengths the hydrodynamics of the effective fluid is fully coupled to 
the elastic dynamics of the worldvolume. For instance this is case for perturbations of thin 
black rings with wavelength comparable to the ring radius. These lie beyond the range of 
applicability of our results. 

It should be quite interesting to extend our analysis to include the extrinsic aspects of the 
blackfold. To do this, one first allows the worldvolume metric where the fluid lives to be a 
curved background, with an extrinsic curvature radius much larger than . This curvature 
acts as an external force on the fluid [3]. In the derivative expansion, the stress tensor 
will in general contain, besides the viscosities, higher-derivative coefficients that multiply 
derivatives of the worldvolume metric. These coefficients will be determined by demanding 
horizon regularity of a perturbation that curves the asymptotic geometry. Perturbations 
of this kind have been studied for certain illustrative examples in [21^ [22| [5] in stationary 
situations that do not involve viscous dissipation. Thus it may be possible to extract the 
extrinsic pressure coefficients in the stress tensor. 

In the AdS context, the external force on the fluid from a worldvolume curvature has 
been studied in [23]. However, in that case the worldvolume geometry is regarded as a fixed, 
non-dynamical background. Instead, in the blackfold context this geometry is dynamical. A 
solution of the forced fluid equations will backreact on the background spacetime where the 
blackfold lives, and thus modify the worldvolume geometry. Therefore for a generic, curved 
blackfold the explicit construction of perturbative metrics becomes rather more complicated 
than in the fluid/ AdS-gravity correspondence. 
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